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We present an information geometric analysis of entanglement generated by an s-wave scatter- 
ing between two Gaussian wave packets. We conjecture that the pre and post-collisional quantum 
dynamical scenarios related to an elastic head-on collision are macroscopic manifestations emerg- 
ing from microscopic statistical structures. We then describe them by uncorrelated and correlated 
Gaussian statistical models, respectively. This allows us to express the entanglement strength in 
terms of scattering potential and incident particle energies. Furthermore, we show how the entan- 
r~{ • glement duration can be related to the scattering potential and incident particle energies. Finally, 

y— -} ' we discuss the connection between entanglement and complexity of motion. 
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PACS numbers: Probability Theory (02.50.Cw), Riemannian Geometry (02.40.Ky), Complexity (89. 70. Eg), 
Entropy (89.70.Cf), Quantum Entanglement (03.65.Ud). 

One of the most important features of composite quantum mechanical systems is their ability to become entangled 
P, In general, quantum entanglement is described by quantum correlations among the distinct subsystems of the 
entire composite quantum system. For such correlated quantum systems, it is not possible to specify the quantum 
state of any subsystem independently of the remaining subsystems |1J. Apart from these remarks, the fundamental 
meaning of quantum entanglement is still a widely debated issue 

From a conceptual point of view, the simplest and most realistic mechanism of generating entanglement between 
two particles is via scattering processes [J, [|| . The two particles can become entangled as they approach each other as 
a consequence of mutual interactions. For instance, for interaction potentials with a strong repulsive core, quantum 
interference between incident and reflected waves can generate transient entanglement. After the collision, the two 
particles may still be entangled and share forms of quantum information in the final scattered state. Quantum 
entanglement can also be generated during inelastic collisions between the dissipative walls of a container and the 
quantum system confined within it [6]. Entanglement may also be induced in multi-atom systems confined in a 
. harmonic trap interacting via a delta interaction potential [7|. 

In order to obtain a clear and detailed understanding of entanglement, it is first necessary to quantify it. It is 
known that for pure maximally entangled states the quantum state of any subsystem is maximally mixed, while for 
separable states it is pure. Thus, one is led to consider the von Neumann entropy of the reduced state, measuring 
its degree of mixedness, as an entanglement measure. This turns out to be correct for pure bipartite states Q (the 
case we are considering in this Letter), while for more general states other entanglement measures should be invoked 
[^. Hoj. Apart from the above presented remarks, a great deal remains unclear about the physical interpretation of 
entanglement measures [lT| and much remains unsatisfactory about our understanding of scattering-induced quantum 
entanglement, especially with regard to how interaction potentials and particle energies control the entanglement 
Finally, our knowledge of the connections between entanglement and complexity of motion remains far from complete 

0. __ 

In this Letter we investigate the potential utility of the Information Geometric Approach to Chaos (IGAC) [13l 1 141 ] 
in characterizing the quantum entangleme nt p roduced by a head-on elastic collision between two Gaussian wave 
packets interacting via a scattering process (l5| . 

IGAC is a theoretical framework developed to study the complexity of informational geodesic flows on curved 
statistical manifolds underlying the probabilistic description of physical, biological or chemical systems. It is the 
information geometric analogue of conventional geometrodynamical approaches flU [ljj where the classical configura- 
tion space is replaced by a statistical manifold with the additional possibility of considering complex dynamics arising 
from non conformally flat metrics (the Jacobi metric is always conformally flat). For recent applications of the IGAC 
to quantum physics we refer to • 

Here we conjecture that the scattering induced quantum entanglement is a macroscopic manifestation emerging from 
specific statistical microstructures. Specifically, using information geometric techniques [20j and inductive inference 
methods [U, , we propose that the pre and post-collisional scenarios are modelled by an uncorrelated [23J and 
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correlated Gaussian statistical model [24| , respectively. We present an analytical connection between the entanglement 
strength - quantified in terms of purity - to the scattering potential and incident particle energies. Furthermore, we 
relate the entanglement duration to the scattering potential and incident particle energies. Finally, we uncover a 
quantitative relation between quantum entanglement and the information geometric complexity of motion (25| . 

Before describing the physical system being studied, we recall that spatially localized Gaussian wave packets are 
especially useful to describe naturally occurring quantum states and they are easy to handle since many important 
properties of these states can often be obtained in an analytic fashion [26j. Furthermore, the Wigner distribution of 
Gaussian wave packets is positive-definite and therefore Gaussian states could be tagged as essentially classical [27]. 

The physical system being considered consists of two interacting spin-0 particles of equal mass m. For such a 
system, a complete set of commuting observables is furnished by the momentum operators of each particle [28j. In 
terms of the center of mass and relative coordinates, the Hamiltonian Ti of the system becomes (Isj . 
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where M = 2m is the total mass and /z = ^ is the reduced mass. The interaction potential V (x) is isotropic and 
has a short range d such that V (x) ps for x > d. Before colliding, the two particles are in the form of disentangled 
Gaussian wave packets, each characterized by a width ero m momentum space. The initial distance between the two 
particles is Rq and their average initial momenta - setting the Planck constant H equal to one - are Tfco, respectively. 
We emphasize that the three-dimensional scattering process in [l5j can be effectively reduced to a one-dimensional 
process as far as the probability density-based analysis concerns. This can be explained by observing that both the 
pre and post-collisional three-dimensional Gaussian wave-packets in (l5j are isotropic. That is to say, in spherical 
coordinates the representations of the separable two-particle states exhibit a functional dependence on the radial 
variable only. For this reason the three-dimensional vectorial representations may be effectively reduced to one- 
dimensional representations. After some algebra, it follows that the initial (pre-collisional) two-particle square wave 
amplitude in momentum space reads, 
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where ±fc are the expected values of the momenta fci and fc 2 , respectively. The square root of the variance for 
each momentum is denoted by oq. Similarly, following [15| and our above mentioned observation on the dimensional 
reduction of the analysis, after some tedious algebra it turns out that the final (post-collisional) two-particle square 
wave amplitude in momentum space in the low energy s-wave scattering approximation becomes, 
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where the adimensional coefficient tqm is defined as, 
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where tqm <C 1 and a s is the s-wave scattering length. As a side remark, we point out that Q reduces to ([3]) when 
^qm vanishes. 

As pointed out earlier, in order to properly analyze entanglement, the entanglement entropy obtained from the 
long time limit post-collisional wave function is required. In most cases however, this must be performed numerically. 
Thus, to approach the problem analytically and simultaneously gain insights into the problem, it is convenient to 
make use of the linearized version of the entropy of the system, i.e. of the purity of the system [l5[. The purity 

function is defined as V = f Tr(p^) where p A = Tr^ {pab) l& the reduced density matrix of particle A and p AB is 
the two-particle density matrix associated with the post-collisional two-particle wave function. For pure two-particle 
states, the smaller the value of V the higher the entanglement. That is, the loss of purity provides an indicator of the 
degree of entanglement. Hence, a disentangled product state corresponds to V — 1. We emphasize that the purity has 
been used as a measure of the degree of entanglement in various physical situations (29| , especially in atomic physics 
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in order to characterize the two-body correlations in dynamical atomic processes |30i. 131], Under the assumption that 
the two particles are well separated both initially (before collision) and finally (after collision) [32| and assuming that 
the colliding Gaussian wave packets are very narrow in the momentum space (ctq <C 1 such that the phase shift can 
be treated as a constant 9 (ko) = 0q), it follows that the purity of the post-collisional two-particle wave function is 
approximately given by [15| 

4 (2k 2 + a 2 ) R y/S^ 
V = l- ° ; + O (S ) , 6 

where So = f 4ir \ f (fco)| 2 is the scattering cross section and / (fco) = e ° fc sm e ° s» jr- is the s-wave scattering amplitude. 

Although very important, the analysis of (l5l | does not address the problem of how the interaction potentials and 
particle energies control the scattering-induced entanglement (a key problem as pointed out in Q) and it does not 
discuss any possible connection between the entanglement generated in the scattering process and the complexity of 
the motion related to the pre and post-collisional quantum dynamical scenarios (a key problem as pointed out in |l2j). 

In this Letter, we attempt to provide some answers to such unsolved relevant issues. We conjecture that the pre and 
post-collisional quantum dynamical scenarios characterized by ^ and ((J), respectively, and describing the quantum 
entanglement (quantified in terms of the purity V in produced by a head-on collision between two Gaussian 
wave packets are macroscopic manifestations emerging from specific underlying microscopic statistical structures. 
We stress that within the IGAC, we provide a probabilist description of physical systems by studying the temporal 
evolution on curved statistical manifolds of probability distributions encoding all the relevant available information 

concerning the system considered. For the problem under investigation, we propose that Pp^^ (fci, fc2|fco, — fco, <jq) 
can be interpreted as a limiting case (initial time limit) arising from an evolving uncorrelated Gaussian probability 

distribution P§£ ) (fci, k 2 \fJ, kl (r) , pi fe2 (r) , a (r)), 

(fci, fc 2 K (r) , M fc2 (r) , o (r)) d ^ != — 2 ±1 (7) 

As a matter of fact, setting /i ki (0) = fco, /U^ (0) = —fco and a (0) = <tq 1 we obtain 

Pg*) (fc 1; fc 2 |fc , - fc , (Jo) = P^f ) (fci, fc 2 |fc , - fc , (To) . (8) 

Similarly, we propose that Pp^P (fci, fc2|fco, — fco, oq\ tqm) can be viewed as a limiting case (final time limit) arising 
from a correlated Gaussian probability distribution Pp^t (fci, fc^l/^ (t) , fJ> k2 (t) , a (r) ; nc), 
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where t"ig = l " klk2 ^ Jf 1 ^^ [ s the correlation coefficient. Indeed, setting fx kl (rg na i) = fco, n k2 (rgnai) — — fco and 
o (rsnai) = cr , we obtain 

pi po2 (fci, felfco, - ^o, co; hg) = fp2f (fci, fc2|fco, - fco, <to; j-q M ) (10) 

in the limit of weak correlations with hq = tqm *C 1. 

At this stage our conjecture is only mathematically sustained by the formal identities (j8]) and ([T0| . To render 
our conjecture also physically motivated, recall that s-wave scattering can also be described in terms of a scattering 
potential V(x) and the scattering phase shift 6 (fc). For the problem under consideration, V(x) equals the constant 
value V for < x < d while it is vanishing for x > d. The quantities V and d denote the height (for V > 0; repulsive 
potential) or depth (for V < 0; attractive potential) and range of the potential, respectively. Within the standard 
quantum framework, integrating the radial part of Schrodinger equation with this potential for the scattered wave 
and imposing the matching condition at x — d for its solution and its first derivative leads to [33] 
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with fc in = y/2fi (T — V) for < x < d, fc out = ^/2^T for x > d and 9q ~ — fco a s denotes the s-wave scattering 
phase shift. The quantities /x and T are the reduced mass and kinetic energy of the two-particle system in the relative 
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coordinates, respectively; k- m and fc out represent the conjugate-coordinate wave vectors inside and outside the potential 
region, respectively. Equation (fTTj) indicates that the scattering potential V(x) shifts the phase of the scattered wave 
at points beyond the scattering region. 

Within the IGAC, given the curved statistical manifolds of probability distributions in and ©, we compute 
the Fisher-Rao information metrics on each manifold. We then integrate the two sets of geodesic equations (three 
equations for each set) leading to the expected trajectories connecting the initial and final macroscopic configurations 
defined by the macroscopic variables f/^, (r) , ^ fea (r) , a (r)). In particular, we obtain that 
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(12) 



Considering Eq. (|T2j) together with the limit of low energy s-wave scattering (kod -C 1) and low correlations (hg <C 1), 
the matching condition (jlip in the information geometric framework reads 



where 



Combining ([13)) and (fT4")l , we obtain 



'(*o) 



riG 



T 



2/jV 



-liVd 3 k . 



(13) 



(14) 



(15) 



Equation (|15|1 . obtained via information geometric dynamical methods, is in perfect agreement with the result pre- 
sented in [34j (and not in [15] ) where standard Schrodinger's quantum dynamics was employed. Such finding is 
especially important because it allows to state that our conjecture is also physically motivated. 

As a consequence of ^ and (|13|) . we find that when both low energy and weak correlation regimes occur, the purity 
V of the system becomes 
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Equation (|T6|) implies that the purity V can be expressed in terms of physical quantities such as the scattering potential 
V (x) and the initial quantities fco, do and Rq via (fl4"|) . Apart from (fT5|) . Eq. is the first significant finding obtained 
within our hybrid approach (quantum dynamical results combined with information geometric modeling techniques) 
that allows to explain how the interaction potential V (x) and the incident particle energies T control the strength 
of the entanglement. The role played by rjc in the quantities V and V suggests that the physical information about 
quantum scattering and therefore about quantum entanglement is encoded in the statistical correlation coefficient, 

specifically in the covariance term Cov (fci, k%) = (fcife) — (^2) appearing in the definition of r\Q. 

Within the IGAC we are also able to estimate the statistical temporal duration over which the entanglement is 
active. Indeed, it turns out that the uncorrelated and the correlated statistical Gaussian models would require time 
intervals r uncorr . and r corr ., respectively, to attain the same value as the initial momentum ko- Assuming 7*1 g <C 1, 
from the above-mentioned geodesic equation analysis it can be shown that what we define "entanglement duration" 
A (fc , er , no ) reads, 
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For example, with 00/^0 ~ 1CP 3 we have Tig < 2 x 10~ 6 . We observe that the entanglement duration can be 
controlled via the initial parameters fc > Co and the correlations riQ (therefore via the incident particle energies and 
the scattering potential due to (fl4|)). Also, we notice that in the absence of correlations, i.e. t\q — > 0, A — > 0. It 
is anticipated that the maximum duration would be obtained when hg is the greatest and the ratio CTo/fco is the 
smallest. In summary, the entanglement duration allows to quantitatively estimate the temporal interval over which 
the entanglement is active in terms of statistical geodesic paths (statistical evolution of probability distributions) on 
the statistical manifolds underlying the information dynamics used to describe the pre and post-collisional scenarios 
(absence and presence of entanglement, respectively). It encodes information about how long it would take for an 
entangled system to overcome the momentum gap (relative to a corresponding non-entangled system) generated by 
the scattering phase shift. The entangled system only attains the full value of momentum (i.e. the momentum value 
as seen in the corresponding non-entangled system) when the scattering phase shift vanishes. For this reason, the 
entanglement duration represents the statistical temporal duration over which the entanglement is active. 

Our final finding uncovers an interesting quantitative connection between quantum entanglement quantified by 
the purity V in (fT6l) and the information geometric complexity of motion on the uncorrelated and correlated curved 



statistical manifolds _/\/f ( uncoir -) anc j A4s CO , respectively. The information geometric complexity (IGC) represents 
the volume of the effective parametric space explored by the system in its evolution between the chosen initial and 
final macrostates. In general, the volume itself is in general given in terms of a multidimensional fold-integral over 
the geodesic paths connecting the initial and final macrostates [25], 
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? 9 n ) is the so-called Fisher density and equals the square root of the determinant of the 



The integration space T>, 



(geodesic) 

e 



V, 



r') in (|2"Tj) is defined as follows, 

ic) (t /) 4ef | Q _ . Q k (Q) < Q k < k (t , } j ) 



(geodesic 

e 



(22) 



(23) 



where k — 1,.., n and 9 



9 k (s) with < s < t' such that, 



d 2 e K ( s ) 

ds 2 



Im 



de l de m 

ds ds 



= 0. 



(24) 



The integration space 

^(geodesic) ^ . n ^ 

is a n-dimensional subspace of the whole (permitted) parameter space 
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The elements of V 
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(t') are the n-dimensional macrovariables {0} whose components 9 are bounded by 
specified limits of integration 9 k (0) and 9 k (t') with k — 1,.., n. The limits of integration are obtained via integration 
of the n-dimensional set of coupled nonlinear second order ordinary differential equations characterizing the geodesic 
equations. In our case n = 3 with = (r) , /i fc2 (t) , a (r)) and following the line of reasoning presented in 
[IE HI one finds that 
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denotes the information geometric complexities of motion on the chosen statistical mani- 
folds. As a side remark, we point out that (|25[) confirms that an increase in the correlational structure of the dynamical 
equations for the statistical variables labelling a macrostate of a system implies a reduction in the complexity of the 
geodesic paths on the underlying curved statistical manifolds |35l.|36|. In other words, making macroscopic predictions 
in the presence of correlations is easier than in their absence. Combining (|T6)) and (f25|) it follows that 
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where. 
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From (|26p it is evident that the scattering-induced quantum entanglement and the information geometric complexity of 
motion are connected. It found that when purity goes to unity, the difference between the correlated and non-correlated 
information geometric complexities approaches zero. In particular, our analysis allows to conclude that the information 
geometric complexity of motion (viewed as an indicator of how "difficult" is to make macroscopic predictions) of 
Gaussian wave-packets decreases when the quantum wave-packets becomes entangled. This is reminiscent of the fact 
that quantum entanglement is generally considered a very useful resource at our disposal in quantum computing. As 
a side remark, we would like to emphasize that the quantification of the complexity of quantum motion represents a 
quite delicate task and it is also a very debatable issue. Indeed, there is no unique manner in which such complexity 
can be described. In particular, the notion of information geometric complexity employed in this Letter differs from 
that introduced in [12| where the number of harmonics of the Wigner function was chosen as a suitable quantum 
signature of complexity of motion. 

In conclusion, our information geometric characterization gives a novel probabilistic picture of quantum entan- 
glement. Within our framework, the key-feature is the study of the statistical evolution of classical probability 
distributions on curved statistical manifolds underlying the information dynamics used to describe the pre and post- 
collisional scenarios (absence and presence of entanglement, respectively). The emergence of entanglement manifests 
itself with a change in the probability distribution that quantifies our state of knowledge about the system. Once this 
is accepted (indeed we provide both mathematical and physical supports to our conjecture, see Eqs. ([8]), ([TO]) and 
(|13p . respectively), our analysis becomes relevant for the following reasons: first, we are able to quantitatively express 
the entanglement strength, quantified by purity, in terms of scattering potential and incident particle energies (see Eq. 
([HJ) ). The scattering potential and incident particle energies are in turn related to the micro-correlation coefficient 
7"ig, a quantity that parameterizes the correlated microscopic degrees of freedom of the system (see Eq. (|14p ). Second, 
we introduce a new quantity termed "entanglement duration" which characterizes the temporal duration over which 
the entanglement is active and show that it can be controlled by the initial momentum pa, momentum spread <To 
and t"ig (see Eq. (|17l) b Finally, we uncover a quantitative relation between quantum entanglement and information 
geometric complexity (see Eq. (|26p b We also point out that our analysis allows us to interpret quantum entanglement 
as a perturbation of statistical space geometry in analogy to the interpretation of gravitation as perturbation of flat 
spacetime. The nature of the perturbation of statistical geometry is global. This, together with the time-independence 
of the geometry, leads to the notion of non-locality. The perturbation of statistical geometry is associated with the 
scattering phase shift in the momentum space. 

We are confident that the present work represents significant progress toward the goal of understanding the rela- 
tionship between statistical microcorrelations and quantum entanglement on the one hand and the effect of microcor- 
relations on the dynamical complexity of informational geodesic flows on the other. It is our hope to build upon the 
techniques employed in this work to ultimately establish a sound information-geometric interpretation of quantum 
entanglement together with its connection to complexity of motion in more general physical scenarios. 
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